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The possibility of realizing bosonic fractional quantum Hall effect in ultra-cold atomic systems 
suggests a new route to producing and manipulating anyons, by introducing auxiliary bosons of a 
different species that capture quasiholes and thus inherit their non-trivial braiding properties. States 
with localized quasiholes at any desired locations can be obtained by annihilating the auxiliary 
bosons at those locations. We explore how this method can be used to generate non-Abelian 
quasiholes of the Moore-Read Pfaffian state for bosons at filling factor v = 1. We show that a 
Hamiltonian with an appropriate three-body interaction can produce two-quasihole states in two 
distinct fusion channels of the topological “qubit.” Characteristics of these states that are related 
to the non-Abelian nature can be probed and verified by a measurement of the effective relative 
angular momentum of the auxiliary bosons, which is directly related to their pair distribution 
function. Moore-Read states of more than two quasiholes can also be produced in a similar fashion. 

We investigate some issues related to the experimental feasibility of this approach, in particular, 
how large the systems should be for a realization of this physics and to what extent this physics 
carries over to systems with the more standard two-body contact interaction. 
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I. INTRODUCTION 

Emergent particles in two-dimensional systems can, in 
principle, transcend the dichotomy of boson and fermion. 
The existence of the exotic statistics stems from the fact 
that winding one particle around another in two dimen¬ 
sions is topologically inequivalent to a process in which 
none of the particles move at all. Thus braiding particles 
does not necessarily bring the system back to the same 
state. Particles obeying fractional braid statisticJ^^ are 
called “anyons.” For these particles, the wave function 
acquires, when one particle winds around another, a 
phase that is a non-integral multiple of 27r. Since the 
product of phase factors is a commutative operation, 
this type of statistics is called Abelian braid statistics, 
associated with one-dimensional representations of the 
braid group. Non-Abelian statistics can arise if there 
is a degenerate set of wave functions for quasiparticles 
at fixed positions. Braiding of some particles around 
the others then corresponds to a unitary matrix trans¬ 
formation in the space of such states. If such matrices 
are non-commutative, the particles are said to exhibit 
non-Abelian statistics. Systems whose excitations sat¬ 
isfy non-Abelian braid statistics are of interest both for 
a demonstration of this physics, and because of their po¬ 
tential application in storage and processing of quantum 
informatiorP. If such degenerate states with quasipar¬ 
ticles at fixed locations are separated from the rest of 
the spectrum by an energy gap, then for sufficiently slow 
perturbations, dynamics is restricted to within such a de¬ 
generate subspace. Any local perturbation has no non¬ 
trivial matrix elements within this degenerate subspace, 
and therefore braiding of quasiparticles is essentially the 
only way to perform nontrivial unitary operations on this 


subspace at low energies, immunizing the system against 
decoherence.l^ 

Of the many theoretical proposals for the realization 
of non-Abelian anyon^SHUl^ the fractional quantum Hall 
(FQH) state of electrons at filling v = 5/2 is perhaps 
the most promising. It is believed to be described by 
the Moore-Read (MR) Pfafian wave functiorP, which 
represents a chiral p-wave paired stat^^^ of composite 

fermion^iSIt^o] 

Numerical studies have shown support 
for this interpretation of the 5/2 FQH state: the exact 
Coulomb ground state at z/ = 1/2 in the second Landau 
level (LL) is in the same universality class as the MR 
state for systems of up to 18 electronsand the MR 
wave function has a significantly lower energy than the 
spin-polarized or the spin-unpolarized composite-fermion 
Fermi sea state in the second LlPl. Abrikosov vortices of 
a chiral p-wave superconductor carry Majorana fermions 
and obey non-Abelian braid statistic^^, and it has been 
proposed that the quasihole excitations of the MR state 
also exhibit non-Abelian statistic^, with 2n quasiholes 
at fixed positions spanning a 2"“^-dimensional degener¬ 
ate spacers. 

Quantum Hall interferometrji^MUl in the solid state 
devices has been the best explored among the propos¬ 
als for demonstration of non-Abelian statistics. It relies 
on the idea that the transport of quasiparticles in the 
edge-currents of a Hall bar is equivalent to braiding them 
around the quasiparticles in the bulk, and the phases ac¬ 
quired in such braiding processes can be detected through 
interference between currents along different paths. Ex¬ 
perimental progress has been made toward the detection 
of such interference^^sHUl^ interpretation of the 

experiments in terms of non-Abelian anyons has not been 
unambiguous. Theoretical work has suggested concep- 
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tual difficulties stemming primarily from the ‘volatility’ 
in the location of the edge and significant effort has gone 
into understanding all the effects that determine the edge 
fluctuations and into extracting signals of non-Abelian 

statistics PStSl 

In this work, we show that some of these problems 
can be overcome, in principle, in cold atom systems. 
Bosonic atoms are neutral, but they behave as charged 
bosons in a magnetic field when subjected to rotation, 
because, formally, rotati on pla ys the same role as a per¬ 
pendicular magnetic fielcPSMZl The filling factor is then 
given by the ratio of boson density to vortex density. At 
sufficiently rapid rotation (i.e. low filling factor), neu¬ 
tral bosons can be driven into the strongly correlated 
regime of the quantum Hall effectP^^^. Another way 
to mimic magnetic field in neutral atomic systems is 
to generat e arti ficial gauge potentials using atom-light 
interactionj ^^ l ^^ l Quantized vortices have been produced 
with various experimental techniques in atomic Bos^^sHsU 
and Fermi gasea^, and extremely high vorticity has been 
achievecpS*^. Progress towards producing the fractional 
quantum Hall effect (FQHE) conditions has been made 
by employing an adiabatic pathway to transfer the atom 
clusters from zero angular momentum in to a g round state 
of any higher total angular momentum L!— I— Bn addition 
to the schemes described above, there have been many 
other proposal^^i®^ aiming to realize FQHE or analogous 
behavior in lattice models, with ultracold atoms confined 
in optical lattic#ill^. In this paper, we focus on contin¬ 
uum atomic gases in two dimensions. While FQHE has 
not been produced in such systems so far, we will as¬ 
sume that it is possible. We will also assume that it is 
possible to engineer arbitrary 2- or 3-body interactions in 
cold atom systems. A method for realization of 3-body 
interaction has been proposed by Roncaglia et re¬ 
cently with implementation requirements within present 
or planned technologies. 

One advantage of cold atoms systems is that it is pos¬ 
sible to introduce auxiliary bosons of a different species 
and manipulate them.ESHSSl’We have shown in a previous 
papei)^ that for appropriate interactions, the auxiliary 
bosons bind quasiholes of the quantum Hall droplet, and 
thus serve as “place holders” for quasiholes. They also 
inherit the fractional statistics of the quasiholes. One can 
thus imagine manipulating the quasiholes by manipulat¬ 
ing the auxiliary bosons. In this article, we construct a 
3-body interaction that accomplishes the same for quasi¬ 
holes obeying non-Abelian braid statistics. We demon¬ 
strate that a specific choice of 3-body interaction can 
be used to prepare states of two quasiholes in either of 
the two different “fusion channels” labeled 1 or ip, which 
differ in their topological character and serve as the el¬ 
ementary qubit of a putative topological quantum com¬ 
puter based on the FQHE. (Swapping the two states can 
be achieved by braiding a third quasihole around one of 
the two quasiholes.) The two states can be distinguished 
through the different Berry phases produced when one 
quasihole encircles another. 


Another advantage of cold atom systems is that sig¬ 
natures of non-Abelian statistics can be seen in the ef¬ 
fective relative angular momentum of a pair of auxiliary 
boson^^. We show that the fractional part of the relative 
angular momentum, which can be deduced from the pair 
correlation function of the pair of auxiliary bosons, corre¬ 
sponds to the fractional Berry phase accumulated when 
one quasihole is taken around another. Using explicit 
numerical calculation of the pair correlation function, we 
show that the fractional angular momenta indeed cor¬ 
rectly reflect the fractional Berry phases associated with 
different fusion channels as predictecP by conformal field 
theory (CFT). This demonstration of non-trivial statis¬ 
tics, while still being at the level of a thought-experiment, 
has the advantage that all the physics happen away from 
the edge of the droplet. Another aspect of non-Abelian 
excitations are the multiplicities of states with excita¬ 
tions at fixed locations. Such degeneracies are reflected 
in the states with auxiliary bosons carrying quasiholes. 

While the situation is relatively clear, at least theoreti¬ 
cally, for a model 3-body interaction, such an interaction 
will not be the easiest to impleme nt in c old atom systems. 
Fortunately, it has been founcP^t^^MI that even bosons 
with 2-body contact interaction exhibit the chiral p-wave 
paired Pfaffian state a.t = 1. Our detailed studies show, 
however, that the applicability of most of the above ideas 
does not extend to a system of bosons with 2-body con¬ 
tact interaction, at least for system sizes available in our 
studies. This is perhaps related to the studied of the 
5/2 state for electrons that show that the excitations of 
the 5/2 state for the Coulomb interaction do not have 
a one-to-one correspondence to the excitations of the 3- 
body interaction for small systems amenable to numerical 
studies. We note that quasiparticles in the MR state have 
also been described theoretically in a way analogous to 
the quasiholes.^sMUl This work focuses on the quasiholes 
of the MR state. 

There have been many other theoretical ap- 
proache^SSKloo] towards the goal of engineering anyons 
in quantum Hall regime and probing their peculiar 
statistics with the help of ultra-cold atoms and optical 
control. Paredes et al^ proposed to create quasiholes 
by focusing laser beams onto the atomic cloud and 
braiding them by adiabatically moving the laser beam. 
Julia-Diaz et a/.l^ and Grab et demonstrated by 
calculating the Berry phase that Abelian fractional 

statistics manifest in ven^ small systemSj _ consistent 

with our previous results^. Kapit et reported 

on a numerical experiment to test the non-Abelian 
braiding properties of lattice bosons, which could be 
implemented with cold atoms in a deep optical lattice 
and an artificial gauge field. Different from the above 
approaches that require time-dependent pot ential to 
braid the quasiparticles. Cooper and SimorP^ proposed 
to demonstrate Haldane’s fractional exclusion statistics 
of quasiholes in the bosonic Laughlin state through 
spectroscopic measurements. 

The organization of the paper is as follows. In Sec. [T^ 
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we review the relevant ideas regarding the MR Pfafiian 
state and its quasihole excitations; and construct wave 
functions describing the 2-quasihole qubit. We then con¬ 
struct a 3-body Hamiltonian in a system of N majority 
bosons and two auxiliary bosons in Sec. Ill which pro¬ 
duce the 2-quasihole qubit states as its exact eigenstates. 
Assuming a system in such a qubit state, we demonstrate 
in Sec. m how the non-Abelian statistics can be seen 
by measuring the pair correlation function of the pair of 
auxiliary bosons. We further study how to produce 2n- 
quasihole states in Sec. [V|and discuss the experimental 
feasibility of our proposals in cold atom systems in Sec. 


VI Finally we conclude in Sec. VII 


II. BACKGROUND 

We begin with a review of the MR Pfafiian ground 
state wave functioiP, the many-quasihole states, as 
well as the explicit qubit representation for 4-quasihole 
state^^. We show that these wave functions can be ex¬ 
pressed in a bipartite form. In particular, two distinct 2- 
quasihole wave functions of the bipartite form are shown 
to be related to the 4-quasihole qubit states with two of 
the quasiholes sent to infinity. 


A. Moore-Read state and quasiholes 


Moore and ReatP constructed a set of trial wave func¬ 
tions for various FQH states using the analogy between 
the correlators in various CFTs and the incompressible 
FQH states. The MR state of the form 


4- 


Pf 

1 / 


m 


Hiz, - z,rpi 

i<j 



( 1 ) 


has a filling fraction ^ where m is even for fermions 
and odd for bosons. Here Zi = Xi + lyt label the com¬ 
plex coordinates of the particles, and £ = ^JhcjeB is 
the magnetic length. This state represents a chiral p- 
wave superconducting state of composite fermionJi^. It 
is also the highest density exact zero-energy eigenstate of 
a 3-body repulsive contact interaction, discussed in the 
following. 

Moore and Read also suggested half-flux quantum ex¬ 
citations of the Pfafhan state described by the trial wave 
function: 

^'2qh OC Y]_{zi - Zj)'^ X Pf(M) X ^ 

i<j 

where M„ = ("■ -- 1^) + (»■ -- hi) 

Zi — Zj 

( 2 ) 

This state has two quasiholes localized at rji = x\ -I- 
ly'i and 772 = -I- * 7/2 > each with a charge ^. (In the 

y- l^.l" 

following we will suppress the Gaussian factor e 
for ease of notation.) 


The wave function for the state consisting of 2n quasi¬ 
holes can be constructed by modifying the argument of 
the Pfafhan in Eq. ([^ as followJ^ 


M ■ = ~ - r]p)...{zj - T]p){zj - + {i^ j) 

Zi Zj 

.. 

In this expression, the 2n quasiholes have been divided 
into two groups of n each (i.e. a, /3 • • • and p, a ■ ■ ■)■ 
There are ways of making such a division. Nayak 

and WilczelP^ demonstrated that they span a space of 
only 2"“^ linearly independent states, consistent with 
the predictions from CFT. This degeneracy can be un¬ 
derstood by noting that each quasihole supports a ma- 
jorana fermion (a half composite fermion). Each pair of 
Majorana fermions in such a state can fuse in two dis¬ 
tinct ways associated with opposite fermion parity, thus 
describing a Hilbert space of dimension 2". Taking into 
account the hxed parity of the overall wave function, this 
produces 2"“^ distinct states. 

A set of linearly independent basis states {pf,a} for 
the space o f 2n -quasihole states has been written by Read 
and RezayP^ as follows: 


•£^F,\{z,rii) = J'^A 



X 


N-F 

2 

n 


^{zf+21-1,Zf+21, V) 
ZF+21-1 — Zf+21 


( 4 ) 


n 

where $( 21 , 22 , 77 ) = X! - 7?r(2i-i)) 

rGS2„ 1=1 


J is the Jastrow factor, A is the antisymmetriza- 
tion operator and S 2 n is the group of permutations of 
{1, 2, 3,..., 2n}. For given locations 
quasiholes, the possible wave functions are indexed by an 
integer F and a sequence A. F is an integer of same par¬ 
ity as the number of particles N, such that 0 < F < n. X 
is an ordered sequence of F integers such that Xi < A^+i 
and 0 < Ai < n. For example, for 2n = 4 quasiholes 
in even N system, indices (F, A) can take two values 
(F = 0, A = [ ]) and (F = 2,A = [0,1]). In general, 
for each choice of F, there are ()() = possible 

choices of A. For 2n quasiholes in even N system, there 
are Z)F=o.even (f) = 2”“^ states for a given {pj. With 
fixed {r]i}, the MR quasihole states shown in Eq. span 
an identical space as that of {(^f.a}- Each MR quasihole 
state can be expressed as a linear combination of the 
basis {:/5 f,a}, in which the coefficients are generally func¬ 
tions of {rji}. These two bases are no longer equivalent 
when {rji} represent positions of auxiliary bosons instead 
of localized quasiholes, as we will consider later, in which 
case only {</5 f,a} are valid wave functions because they 
are symmetric in the Ty^’s. (The MR 2-quasihole state 
is an exception because it is already symmetric under 
exchange of rji and 772 -) 

The above discussion gives the counting for quasihole 
states in a system without any edge degrees of freedom or 
equivalently in a spherical geometry. In the presence of 
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an ‘open’ edge, the parity constraint is relaxed because 
the edge can accommodate a majorana fermion. This 
is in addition to the purely bosonic excitations which 
are fully localized at the edge. In the above mentioned 
basis <^F,A) the upper limits on values of F and A arise 
from a hard angular momentum cut off representing the 
edge. These constraints are again relaxed in the presence 
of edge degrees of freedom, thus resulting in increased 
degeneracy. 

Multiple degenerate states with fixed quasihole posi¬ 
tions are necessary for non-Abelian statistics. The sim¬ 
plest case with such a degeneracy has four quasiholes, 
whose states form a two-dimensional Hilbert space.l^ A 
basis set can be chosen following Nayak and 

Wilczek which has simple braiding properties: 





(?? 13 ? 724 )^ 

(l±V^)l /2 


( 5 ) 


X ('I'( 13 )( 24 ) ± V^'f'( 14 )( 23 ))e ■ 

where x = ryi4??23/?7i3?y24 and 77^/3 = rja - rip. In this 
expression. 


i<j 

where M ~ ~ ~ ~ ^ j) 

Zi Zj 

( 6 ) 

An adiabatic transport of the quasihole 771 around 772, 
or of 773 around 774 results in an overall Abelian phase 
depending on the state (I or tjj). In contrast, braiding 
of either 771 or 772 around either 773 or 774 results in an 
interchange of the states and 4 '’^ as a result of the 
branch cuts in the ^/x factor. If the two states represent 
a qubit, then such braiding operations serve as simple 
gates. 

The states and 4 '’^ correspond to conformal blocks 
of a correlator in an Ising CFT as described later in this 
paper. The quasihole insertion operators correspond to 
the cr fields and the two states above correspond to the 
fusion of cr( 77 i) and <7(772) into 1 or 7/1 channels. 

Based on a plasma analogy, Bonderson, Gurarie and 
Nayak demonstrated that and vjiCb) ^re orthogonal 
in the lim it where the quasiholes are far apart from each 
other.ESH Here we directly compute the overlap between 
the states 47(^1 and by Monte-Carlo technique for 
many different quasihole positions (chosen randomly) for 
a relatively large system N = 60 . Note that all the over¬ 
laps calculated in this p aper are norm alized, denoted 
as {ipi\\ipj) = Figure 0 (a) 

shows the amplitude of the overlap ( 47 *^^^ 1147('‘^)) as a func¬ 
tion of the minimum distance between two quasiholes in 
each set of quasihole positions. The calculation shows 
that the overlap is zero (within numerical uncertainty) 
when the quasiholes are far enough from each other (e.g. 
\''li~Vj \ ^ 2 ). When two of the quasiholes approach each 
other, there are two branches shown in the figure, where 
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FIG. 1: (Color online) (a) The modulus of the overlap be¬ 
tween the 4-quasihole qubit states 47 ( 7 ) vj/l’/’) g^g g func¬ 

tion of the minimum distance between two quasiholes in each 
quasihole-position configuration. The blue squares represent 
cases when the minimum distance is either |r 7 i — 772 ! or I 773 — 774 I, 
and red squares are other configurations, (b) The modulus of 
the overlap between the 2 -quasihole qubit states 47 (i)( 2 ) and 
47(1)(2) — 247 (i 2){) as a function of the maximum quasihole 
position |77i|/i?jv where Rn ~ \f2Njv is the radius of the 
quantum Hall droplet. In both calculations the quasihole po¬ 
sitions are chosen randomly, and the system size is A = 60. 


the red squares represent the configurations in which ei¬ 
ther 771 or 772 is close to either 773 or 774, and the blue 
ones indicate 771 is close to 772 or 773 is close to 774. We 
thus conclude that the overlap between 47 ^^) and 47 (’(') is 
essentially zero except when 771 or 772 approaches 773 or 
774. 


B. Bipartite representation and 2-quasihole qubit 


The Pfa ffian wa ve function can be rewritten in a bi¬ 
partite forn l^°^b ^ 


Ar/2 Af/2 

'I'yil/m = J^(z2j - Zij)^ (2:23-1 - 2:21-1)^, 

N 

j= “ ^ 7 ) 

i<j 


( 7 ) 
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where S is the symmetrizer over all the complex coor¬ 
dinates {zi}i=i^N■ The relation to the Pfaffian wave 
function can b e derived with the help of an identity due 
to CauchjIiSl. 

Our study focuses on the bosonic Pfaffian states at 
filling V = \ (namely m = 1 ). This wave function is ob¬ 
tained by symmetrizing two partitions of incompressible 
states at filling fraction v = Quasihole wave functions 
can be produced by now inserting a single quasihole into 
each of these partitions as follows 


^(l)(2) - 




N/2 N/2 

where ;r?) = (cci - 77) - xj)'^ 

i=l i<j 

( 8 ) 

The subscript (•••)(■••) of dt on the left hand side, in¬ 
dicates the quasihole coordinates in the two partitions. 
Specifically, ( 1 )( 2 ) implies that r]i and 772 belong to differ¬ 
ent partitions, as opposed to (12)() that refers (described 
below) to a state in which both quasiholes reside in the 
same partition. We demonstrate in Appendix that 
this wave function is exactly the same as the MR 

2 -quasihole wave function written in a Pfaffian form as 
shown in Eq. ([^. The bipartite representation immedi¬ 
ately suggests another 2-quasihole state in which the two 
quasiholes are placed in the same partition: 


N/2 N/2 

^(12)0 = S W{z2i - m)(.z:2i - 772) J|(2:2i - Z2jf 

i —1 i<.j 

N/2 

X (2:27-1 ~ 22^-1)^. 

i<j 


( 9 ) 


This corresponds to the Pfaffian wave function for the 
2 -quasihole state (see Appendixfor details): 

W(12)()=Pf| -- [[[Z^-ZJ} . 

\ '^3 / 

( 10 ) 

The two different 2 -quasihole states 'I'(i)(2) and 'I7(i2)() 
are related to the two different 4 -quasihole states in Eq. 
([^. More precisely, the space spanned by these two 2 - 
quasihole states can be obtained by taking two of the 
four quasiholes in Eq. © to infinity by setting 773 = r, 
774 = re*® and taking the limit r —>■ 00. In this limit, we 
find (see Appendix [B|) 

^ = (ffi - ^2)'/®'k(l)(2) 

VI/W^ vI/W=(r7i-772)^/8(vI/(i)(2)-247(12)0). 


Therefore 4'(i)(2) and 4'(i)(2) — 247 (i 2 )o form a 2 -quasihole 
representation for the qubit states. Note that these re¬ 
duced states do not transform into one another under any 
braiding operation; braiding 771 around 772 simply pro¬ 
duces an Abelian phase. Nevertheless, the fusion prop¬ 
erties of the 771 and 772 quasiholes remain the same as 


those in the 4 -quasihole qubit system: 771 and 772 fuse into 
identity 1 in 4'(i)(2), while in 4'(i)(2) — 24 '(i 2 )() they fuse 
into The fusion properties are reflected in the differ¬ 
ent Berry phases associated with braiding one quasihole 
around another, which we discuss in detail in Sec. IV 
The 2 -quasihole states thus form a two-level system, i.e. 
a qubit. The state can be flipped by introducing a third 
quasihole that braids around one of the two quasiholes. 
This is exactly the picture of the constricted Hall bar as 
a quantum bit for making fractional quantum Hall quan- 
.|tum computer.l^^t^ 

In a finite system, the two quasiholes that are taken 
to large distances actually reside at the edge of the sam¬ 
ple and thus the two wave functions differ at the edge of 
the system. This can also be seen when the above wave 
functions are realized in the spherical geometry. When 
placed on a sphere with a flux of 2 Q = m{N — I), the 
wave functions 'I'(4)(2) (47(12)0) bas two additional quasi¬ 
holes in opposite partitions (same partitions) at the south 
pole. Additional edge excitations are possible but do not 
influence the properties of the bulk quasiholes. 

As discussed above, the 4 -quasihole qubit states are 
orthogonal for well-separated quasihole positions. The 
same is true for the 2-quasihole qubit states 47(i)(2) 
and 47 (i((2) — 247 (i 2 )(). For relatively large system size 
[N > 60 ), we have performed many trial computations 
for different sets of rji positions, and found the overlaps 
are essentially zero except when a quasihole is located 
close to the edge of the FQH droplet. The results are 
shown in Fign(b). Considering that 7734 actually reside 
at the edge, the results indicate that the overlap is non¬ 
zero when either of 771 and 772 is close to 773 or 774, and is 
zero otherwise, independent of the distance It/i— 772]. This 
is consistent with the results for the 4 -quasihole qubit 
states shown in Fig. 0(a). 


III. PRODUCING THE QUBIT STATES WITH 
3-BODY INTERACTION 

A. Quasiholes versus auxiliary bosons 

The previous section discussed the nature of localized 
quasihole excitations of the Pfaffian state. In describing 
the quasiholes, the locations are mere parameters. In our 
earlier worlPSl^ we proposed to introduce a pair of aux¬ 
iliary bosons into a background FQHE state of bosonic 
atoms. The auxiliary bosons capture localized vortex¬ 
like quasihole excitations of the quantum Hall state and 
then behave as particles with fractional braid statistics. 
In fact, a system with two species of particles is a nat¬ 
ural platform for producing quasiholes, where the parti¬ 
cles of one species (namely the auxiliary bosons) serve as 
“place holders” for the quasiholes of the other species. In 
other words, the quasiholes of the FQH state are enslaved 
by the auxiliary bosons. This idea can be illustrated 
with wave functions. For example, the 2 -quasihole state 
47(4)(2) in Eq. ^ turns into a wave function describing 
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two species of bosons by changing the quasihole-position 
parameters ry^’s into dynamical coordinates and including 
a Gaussian factor for them: 

Njl Njl N/2 

'I'(l)(2) = S - 111) ” %) 

i=l i<j i-1 

NI 2 

i<j 

( 12 ) 

where S symmetrizes the z-bosons. Here the ryi-bosons 
and the z^-bosons share the same magnetic length as re¬ 
flected in the Gaussian factor. Note that 'I'(i)(2) is al¬ 
ready symmetric under exchange of rji and 772, making 
it a valid wave function even when the Ty^’s are thought 
of as auxiliary boson coordinates. That is also true for 
\l'(i2)(), but not for MR quasihole states with 2 n > 2 . 
In the following, we will use and 'I'(i2)() to de¬ 

note the corresponding wave function for two species of 
bosons, which is slightly different from the wave function 
for localized quasihole states in Eqs. Q and ([^, but the 
meaning should be clear from the context. 

We ask if there exists a Hamiltonian that produces the 
qubit states 'I'(i)(2) and 'I'(i)(2) ~ 2'I'(i2)() as its exact 
eigenstates. Before proceeding to that question, we first 
study the orthogonality and normalization properties of 
the 2 -quasihole states. 


B. Orthogonality of 2-quasihole qubit states 

We have written down two different 2 -quasihole states 
with the help of bipartite representation in the last sec¬ 
tion as 'I'(i)(2) and 'I'(i2)()- In order to study their prop¬ 
erties in the thermodynamic limit, we calculate the ratio 
of the normalization factor of the wave functions 'I'(]^)(2) 
and ^'(12)0 




N N 



FIG. 2: (Color online) (a) Black square represents the ratio 
of the normalization factor of the two different 2-quasihole 
wave functions as shown in Eq. (131 and red circle represents 

(HI 


the normalized overlap between them as shown in Eq. 
as a function of system size N. (b) The natural log of the 
corresponding quantities as a function of N. In panels (a) 
and (b), we use ||i/>i|P to denote the normalization constant 
{ipilipi). (c) The normalized overlap between the 2-quasihole 
qubit states (5'(i)(2)| |^'(i)( 2 ) - 2’!' (i 2 ){)) US a function of N. In 
all panels ryds have been taken as dynamical coordinates in 
Monte-Carlo calculation. 


We also compute the normalized overlap between the 
two 2 -quasihole qubit states: 


(^(l)(2)l^(l)(2)) 
(^( 12 )()|^( 12 )()) ’ 

as well as the normalized overlap 
(^(l)(2)l^( 12 )()) 

\/ ('I'(l)(2)l'I'(l)(2))('I'(12)()l'I'(12)()) 


( 13 ) 


( 14 ) 


for different system sizes. Note that we have taken the 
ryds as the dynamical coordinates of the auxiliary bosons 
in all the calculations. Fig. (a) shows that the values 
of the two quantities decrease as the number of particles 
N increases. The plot of the natural log of the quantity 
versus N in Figure [^b) indicates that both Fqs. ( 13 ) 
and ( 14 ) decrease exponentially with N. In the limit 
N ^ 00, 4 '(i 2 )() dominates in the expression 4 '(i)( 2 ) — 
2 'I'(i 2 )() and is orthogonal to 'I'(i)(2). We thus conjecture 
that in the thermodynamic limit, the second qubit state 
4 '(i)( 2 ) - 24 '(i 2 )() reduces to 4'(i2)()- 


(4/(i)(2)|4'(i)( 2) - 24f(i2)()) 

V(^(1)(2)I'I'(1)(2))*('I'(1)(2) — 24'(i2)()|^(1)(2) “ 24'(i2)()) 

( 15 ) 

Figurel^c) shows that the overlap is very small (< 0 . 07 ) 
even in small system {N = 6 ) and becomes essentially 
zero for N > 30 . This result implies that 4 '(i)( 2 ) and 
'k(i)(2) — 2 'I'(i 2)() serve as a good orthogonal basis for the 
2 -fold degenerate space of two quasiholes. This result 
will be useful in the next subsection where we show that 
both 'k(i)(2) and 4 '(i)( 2 ) — 2'I'(i2)() are exact zero-energy 
eigenstates of a certain 3 -body Hamiltonian. 


C. Model Hamiltonians 


The Pfafhan ground state is known to be the ex¬ 
act zero-energy state of the repulsive 3 -body interaction 
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FIG. 3: (a) Energy spectrum obtained by diagonalizing the 
interaction V = I4zz + I4zrj within a Hilbert space formed by 
N = 8 jz-bosons and two 77 -bosons with L = 32 and = 

12. There are two zero-energy states in the spectrum, (b) 
Energy spectrum obtained by adding the 0.01I4>7>7 term to 
the interaction in (a) and diagonalizing the interaction within 
the same basis. 

Hamiltoniaiii 22 I 

N 

H = Vo 5 {ri - rj) 5 {ri - fk) ( 16 ) 

i<j<k 

We consider a system of two species of bosons: N z- 
bosons and two 77-bosons, and take the 3 -body interaction 
to be identical for all particles. In this sense, we write the 
interaction as y = Vzzz + Vzzrj + VzTjr]- We perform exact 
diagonalization within the basis of total angular momen¬ 
tum L = 7 V^/ 2 , corresponding to the 2 -quasihole states. 
Note that and ^'(12)0 have the same total angu¬ 

lar momentum L, but different maximum single-particle 
angular momentum for the z-particle: Zmax = N — 1 and 
N respectively. Here and in the following we consider 
a, N = 8 system and take the cut-off on single-particle 
angular momentum /max to be 12 for both z and 77 par¬ 
ticles. (In fact we have tested that, as long as /max > N, 
the value of /max does not affect the zero-energy states 
obtained by diagonalizing any of the 3 -body Hamilto¬ 
nian we have considered here and following in the 2- 
quasihole case.) The exact energy spectrum shows only 
one zero-energy state under this interaction, which is ex¬ 
actly «'(l)(2)- 

In order to produce 'I'(-i2)(), we observe from the bipar¬ 
tite representation that the interaction V = Vzzz + Vzzr) 
annihilates the wave function of d7(]^2)() as well as d7(]^)(2)- 
We thus diagonalize V = Vzzz + Kzr; in the same basis as 
above and obtain two zero-energy states in the spectrum, 
shown in Fig. (a). The 2 -fold zero-energy subspace is 
spanned by 'I'(i2)() and 'I'(i)(2) because either of the zero- 
energy states can be expressed as a linear combination of 
'I'(i2)() and 'I'(-i)(2)- We then use a small 14 ,,,^ interaction 
to lift the degeneracy. Particularly, we diagonalize the 


Hamiltonian 

V = Vzzz + Vzzr, + O.OlVzrjr,, ( 17 ) 

and obtain a low-energy doublet including one exact zero- 
energy state and one close-to-zero energy state as shown 
in Fig. 1 ^ (b). The zero-energy state is still exactly 
'i'(i)(2), and the close-to-zero energy state approximately 
corresponds to the orthogonal state in the 2 -fold zero- 
energy subspace. This approximation should become ex¬ 
act when the strength of 14,,,, interaction approaches in¬ 
finitesimally small value. Here for N = 8 , 4 '(i 2 )() is not 
orthogonal to 'I'(i)(2) though they become orthogonal at 
thermodynamic limit as discussed in Sec. m On the 
other hand, the linear combination d 7 (,^ 7 ( 2 ) ~ 2'I'(22)() is 
essentially orthogonal to 'I'(i)(2) even in relatively small 
systems as we have shown in Fig. |^(c). Hence the upper 
state of the low-energy doublet is essentially the second 
2 -quasihole qubit state 'I'(i)(2) — 2 'I'(i 2 )(). In the limit of 
infinitely small 14 ,,,, and large N, the first excited state 
becomes exactly d>(i)(2) — 2 d>(i 2 )(). We thus demonstrate 
that the qubit, namely the two-level system formed by 
two quasiholes, can be produced with the 3 -body Hamil¬ 
tonian V = I 4 zz+I 4 z,,, and each single qubit state can be 
precisely obtained by including a small Vzrjr, interaction 
that lifts the degeneracy between the two states. 

IV. FRACTIONAL BRAID STATISTICS 
A. CFT representation of states 

We briefly outline the derivation of the two qubit states 
for four quasiholes and their braiding properties. By gen¬ 
eralizing the analogy between wave functions for simple 
FQHE states and chiral correlators in certain confor¬ 
mal field theories, Moore and Read proposed the Pfaf- 
fian wave function, which can be written as a correlator 
within an Ising CFT consisting of three primary fields 
1,7/), cr of conformal dimensions 0 , j, ^ respectively to¬ 
gether with a bosonic field cj)W The real fermion/boson 
forming an FQHE state is identified as which is 

fermionic if m is even and bosonic if m is odd. The Pfaf- 
fian wave function made up of N such fermions/bosons 
can be written as the the correlator of N such opera¬ 
tors. The quasihole operator is chosen to be 
which satisfies the requirement that the wave function be 
single-valued in particle coordinates and have the small¬ 
est charge {Q = l/ 2 m). A wave function with localized 
quasiholes is obtained by inserting the quasihole opera¬ 
tors in the correlator. In particular, the wave function 
for a 4 -quasihole state of a system with (even) N fermion 
or boson can be written as 

^'4qh =(ct(?7i) • • • ■ ■ ■ ■fpizN) 

X 

X g*V™</'(zi) . . . ^iVmij>{zN) ^-i J iV Zy/mpo<p{z)\^ 

( 18 ) 






The above correlator decomposes into two conformal 
blocks of {aaacrxj} ■ ■ ■ ip) corresponding to the two pos¬ 
sible ways that the four tr operators can fus e to give 1 
under the fusion rules of the Ising CFT 1^221 The fields 
<j{r]i)a{ri2) and the fields a{ri^)a{r]i) can each fuse to a 
1 or ^ fields and the two choices result in the two wave 
functions and Exact form of t he con formal blocks 
can be obtained through bosonizationPSESU The explicit 
wave function given in Eq. ([^ can be obtained with the 
help of the identitjff^ 


the majority bosons forming the quantum Hall droplet. 
It was shown that the fractional angular momentum re¬ 
sults in a “quantization” of the separation R between the 
auxiliary bosons given by the expression R? ^ 2 Mes£ ^ 
where Meff = M — v and £ ^ , which may be ob¬ 

servable in measurements of the pair correlation function 
for the auxiliary bosons. 

The combined wave function of the two auxiliary 
bosons and the background quantum Hall droplet takes 
the form 


= Ylizi - ZjY 


(19) 


i<j 


Charge neutrality required in the above equation is en¬ 
sured by the smeared background charge.l^ It has been 
conjecturecP^that if the wave functions of states are writ¬ 
ten in the basis specified by the conformal blocks, the 
Berry phase vanishes - all braiding properties are man¬ 
ifested in the basis functi ons. T his can be understood 
using the plasma analogjff^Hinil and has bee n numeri¬ 
cally verified using Monte-Carlo method 4 i 22 l. Thus we 
can now read off the braiding properties directly from 
the basis wave functions. As we have discussed in Sec. 


2 -quasihole states [shown in Eq. 0 ] if two of the quasi¬ 
holes are sent to infinity. Since the Berry phase matrix 
vanishes in this basis, the phase gained when braiding rji 
around 772 should be equal to 


HB 


the 4 -quasihole qubit states reduce to the two 


27 r 

Y 

for 

^(l)(2) 

(20) 

IOtt 

“8“ 

for 

4 '(i)( 2 ) - 24 '(i 2 )() 

(21) 


Here, 47 ( 1 )(2) and 4 '(r 2 )() describe the quasihole states 
rather than the states with two species of bosons. How¬ 
ever, since bosons do not produce any additional Berry 
phases (modulo 27 r) for exchanges or windings, the bound 
state of a boson and a quasihole should obey the same 
fractional statistics as the quasihole. We verify these re¬ 
sults numerically in the following subsections. 


B. Effective relative angular momenta 

The possibility of fractional relative angular momen¬ 
tum was introduced in a previous worli^ in the context 
of quantum Hall states at filling fractions ly = , which 

we briefly review here. The relative angular momentum 
of a pair of bosonic particles in the lowest Landau level 
(LLL) is quantized to even numbers. It was shown that 
if the same bosons bind anyonic quasiholes of a quantum 
Hall droplet, the boson-anyon bound states exhibit frac¬ 
tional effective relative angular momenta whose values 
are directly related to the Abelian fractional statistics 
of the anyons. It was demonstrated that the binding of 
a boson to an anyonic excitation can be achieved with 
a contact interaction between the auxiliary bosons and 


Y{z,m,m) ^ 


N 

'i’iy{z) X ]^(^* - ? 72 )(-Zi - 
2=1 


xe 




IP 


Vi) X Pm{vi,V 2 ) 

( 22 ) 


where '^u{z) is the wave function of the homogeneous 
incompressible quantum Hall state at filling v and 


■Pm(»7i,»72) = XI ‘^L{r]i-V2)^iRi+ri2)^ ^ (23) 

L=even 


is a symmetric polynomial with total angular momentum 
M. In our evaluations of the pair correlation function we 
will assume the most symmetric situation in which 772 lies 
at the origin. Then, with 77 = 771 denoting the relative 
coordinate, the wave function reduces to 


N 

1p{z,T]) ~ 'i>u{z)Y[[z^{Zi 
2=1 


v)]v 




hr 


(24) 


where M can take all integer values. 

The relative distance between the bosons in such a 
wave-function can be estimated using a simple semi- 
classical picture in which the bosons can be thought of as 
revolving around each other at a distance R. The most 
probable value of R is fixed by the requirement that the 
total Aharanov-Bohm phase accumulated along an orbit 
of radius R be equal to the Berry phase gained by the 
wave-function upon winding a boson around the other at 
the same radius. The latter corresponds to the number of 
vortices that rji sees on the other bosons enclosed by the 
orbit. Using a mean field assumption of uniform distri¬ 
bution of bosons in the quantum Hall droplet away from 
the quasiholes, we can arrive at the above mentioned de¬ 
pendence of the radius R. 

The above argument can be rephrased in terms of a 
conjecture that the quantum Hall quasihole wave func¬ 
tions when written as the conformal blocks of appropri¬ 
ate CFTs have a normalization consta nt t hat is indepen¬ 
dent of the locations of the quasihole^^^Sl. For example, 
consider the 2-quasihole wave function at filling fraction 
i' = 1 /q of bosons: 


N 

V'CFT,i(^>? 7 l:'^ 2 ) OC (771 - 772)’ nn (Zi - r]j) 

i=lj = l ,2 


N 

X Yl {Zi 
i<j = l 




E n 
2 = 1 


UiT _ \vir+\r]2r 
4 qP 


( 25 ) 
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For sufficiently separated quasihole locations, 
the conjecture is that the normalization Af = 
f 1 (z, 771 , ry 2 )P is weakly dependent on 7712 - In 

terms of this, the pair distribution function ( 7 ( 771 , 772 ) of 
the combined wave function of the two species of bosons 
in Eq. is given by 


G{r]i,V2)(X- dz\ip{z,r]i,T]2)\" 


_n_i'I 

xAfe ^ ^ . 


(26) 


The peak of the pair correlation function G{r]i — 0, [ 772 1 = 
R) is then located at 


= 2{M - u) 


1 - V 


(27) 


This derivation of the relation between M and R? 
can be extended to the case of the 2-quasihole wave 
functions of the Pfaffian. In the above example of the 
Abelian quantum Hall state, the 2-quasihole state cre¬ 
ated in the quantum Hall droplet due to the repulsion 
between the two species of bosons corresponded to a 
unique wave function with two single-quantum vortex ex¬ 
citations. The Pfaffian presents the possibility of a qual¬ 
itatively different behavior. For a suitably chosen inter¬ 
action between the species (as discussed in the previous 
section), the Pfaffian quantum Hall wave function pro¬ 
duces half-quantum vortex excitations at the locations 
of the auxiliary bosons. Unlike the Abelian case, here 
there are two degenerate states with such quasiholes. 


For the MR quasiholes, the wave function in Eq. (22) 
is replaced by 

'>P{z,Vl,V2) ^ '^2cih{z,Vl,V2) X PMijli.m) (28) 

where 'I' 2 qh(-Zj 7 ii ?? 2 ) could be an arbitrary state in the 
two-dimensional space of 2 -quasihole states 

^ 2 qh(^, 72 ) = + &('I'( 1 )( 2 ) ~ 2\I'(i2)())- (29) 

The only dependence of d 72 qh(z, 771 , 772 ) on rjiS is through 
the wave functions 'k(i)( 2 ),(i 2 )() and not through coef¬ 
ficients a and b. In te rms of normalizations of 

(shown in Eq. (B7|), which by conjecture are 
weakly dependent on 771 2 , we can write the pair correla¬ 
tion ( 7 ( 771 , 772 ) in the wave function in Eq. as: 


G{vi,m) oc [\a\'^\PM{vi,V2)\'^\m - m 
\b\^\PM{m,m)\‘^\m -m 




+ 




The cross term between 'k^qh and 4 ' 2 qh is proportional to 
the overlap between the two states and has been assumed 
to be zero in arriving at the above result. Maximization 
of the pair correlation function (7(0, R) with respect to 
R^ as before gives 


4 M- - 


< < 4 ( M - - 

^ ^ 8 


(31) 


where the right and left hand limits correspond to the 
(a, b) = (0,1) and (1, 0) respectively. The fractional part 
of the effective relative angular momentum, as can be 
read out in Eq. pll), exactly gives the Berry phases 


(l/8)27r and (5/8)27r in Eqs. 20 and 21 
The relation in Eq. (31) can also be derived within 


a simple semiclassical approximation mentioned above 
for quantum Hall states with Abelian anyons. Eor 2- 
quasihole Pfaffian states, one needs to consider an ad¬ 
ditional statistical phase that 771 sees on 772 as¬ 

sociated with different fusion channels. We assume 772 
is fixed at the origin, and equate the number of flux 
quanta {ipo = hc/e) enclosed by the circle of radius i?, 
i.e. 'KR^B/(j)Q = R? , to the mean number of enclosed 
vortices \N M -\- . Here N = vR?I2P — is the 

average number of enclosed z particles, and the | factor 
in front of N comes from the fact that the quasihole of 
Pfaffian state has only half-quantum vortex; the last term 
— reflects the expulsion effect of the centrally located 
quasihole. Solving for i?^, we obtain 


R^ =4(M-h<5(i’’^) -77/4)£2 


(32) 


with V = 1 for bosonic Pfaffian state. The phases 
cannot be estimated from such a semi-classical picture 
but can be seen to be = 1/8 and = —3/8 by 
comparison with the limiting cases of Eq. (31). 

The mean field derivation of Eq. (321 above has as¬ 


sumed the condition that the 77 -particle is inside the disk 
of quantum Hall state of z. For larger M, 771 is pushed 
outside the quantum Hall droplet, and the number of en¬ 
closed vortices is simply M + N/2 + which leads 

to the relation 

i?2 = 2(M -h N/2 + (outside). (33) 


C. Numerical studies 


We next verify the two limiting cases in Eq. (31) 


with explicit calculation of the pair correlation functions 
5 ( 1771 1 = r, 772 = 0) for states specified by (a, 5) = (0,1) 
and (1,0) respectively. We fix one auxiliary boson at the 
origin, and numerically calculate the density profile of 
the other one with Monte-Carlo techniques. As a refer¬ 
ence, we first consider a system containing only the two 
auxiliary bosons 771 and 772 for 2 < M < 14 (with no z 
bosons), in which case the peak of the pair correlation 
function is simply given by R^ = 2M£^^ The result is 
shown in Eig. (a). As we introduce the two auxiliary 
bosons into a bosonic Pfaffian state of z-bosons, the den¬ 
sity profile of the bosons screened by the quantum Hall 
fluid “moves outwards” as shown in Pig. |^(b), where the 
solid line displays the pair correlation for 'l'(i)( 2 ) state 
and the dashed line for (4'(i)(2) ~ 247 ( 12 ) 0 ) state. The 
change in the length scale governing the auxiliary pair 
of bosons (£^ — 2 £^) is mainly due to the reduced ef¬ 
fective magnetic held caused by the interaction between 
auxiliary bosons and background bosons.l^ 
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FIG. 4: (Color online). The pair correlation function g{r,0) 
of the arrxiliary bosons with one boson fixed at the center (a) 
without any background z particles and (b) in the background 
of a = 1 bosonic Pfaffian state of = 20 jz particles. The 
solid and dashed lines in (b) show the pair correlation func¬ 
tion for the different 2 -quasihole states 4 '(i)( 2 )T’Ar(? 7 i,r? 2 ) and 
('l'(i)( 2 ) - 2 'I>(i 2 )())Pm(? 7 i,» 72 ) respectively. The pair correla¬ 
tion functions here and in the following are all quoted in units 
of Po ~ (27r£^)“^, where £ is the magnetic length. 


To further explore the fractional statistics from the 
pair correlation function, we plot the square of the peak 
positions as a function of M and N in Figs, [^a) an d^b) 
for the two states and ~ 2 'km)()). The 

predicted behaviors from Eqs. ( | 32 [ ) and ( | 33 [ ) are also 
shown with black and red dashed lines for comparison. 
The behavior predicted by Eq. ( 32 ) is fully confirmed for 
quasiholes in the interior of the quantum Hall droplet. 
The behavior changes from Eq. ( 32 ) to Eq. ( 33 ) for large 
M, when the auxiliary boson rji lies outside the droplet. 
Note that the fractional part of the effective angular mo¬ 
mentum 7 = — 1//4 is just given by the x-intercept 

of the blac k da shed line. To study 7 more precisely, we 
write Eq. ( 32 ) as B?!£? = 4 (M -|- 7 ) and solve for 7 : 
7 = — AM). We therefore plot the quantity 

\{B? /P — AM) as a function of M for = 40 and 50 in 
Fig. (c). This quantity is displayed for the three wave 
functions: ^'(i)(2), ^(12)() and 4'(i)(2) - 24 '(i 2 )(), and the 
expected values, 7 = —1/8 for 'I'(i)(2) and 7 = —5/8 
for d>(iU2) — 24 '(i 2 )(), are also plotted with dashed line. 
Figure iMc) shows that the numerical result of 7 is con¬ 
sistent with the expected value, and the state 'k(i2)() 
(corresponding to (a, 5 ) = (1,-1)) produces a value 
in between the two limits in Eq. ( 31 ). The deviation 
at small M {M = 1 , 2 ) is unsurprising because quasi¬ 
holes need to be well separated to exhibit the expected 
braid statisticJ^^^HIII!^ and the deviation at larger M for 
iV = 40 indicates proximity of 771 to the edge. Figure 
[^d) plots the differences of the quantities in (c) between 
^(l)(2) and ^'(12)0 (black square), and between 'I'(i)(2) 



FIG. 5: (Color online). The peak positions of the pair cor¬ 
relation function as a function of the relative angular mo¬ 
mentum M for N = 6 — 60, for (a) state and (b) 

^(l)( 2 ) — 2 T(i 2 ){) state. The black dashed line plots the pre¬ 
dicted behavior in Eq. ( |32[ ) and the red dashed line plots Eq. 

f . (c) The quantity j(i?^/l^ — AM) as a function of M 
N = AO and N = 50 for (black square), 4 '(i 2 )() 

(blue circle), and T(i)( 2 ) — 24'(i2)() (red diamond). The black 
dashed lines plot the expected 7 value: —1/8 and 5/8. (d) The 
differences of the quantities |(i?^/Z^ — AM) between 4 '(i)( 2 ) 
and T(i 2 )() are shown in the black square, and between 4 '(i)( 2 ) 
and 4 '(i)( 2 ) — 2 T(i 2 )() are shown with the red circle. 


and 'I'(i)(2) — 24 '(i 2)() (red circle). The latter one is ex¬ 
pected to be 0 . 5 . 

From these numerical results, we stress that even 
though it is not easy to detect the accurate fractional 
Berry phase from small system {N = 6 • • • 10 ), the pair 
correlation of the auxiliary pair has an obviously different 
profile for the two topologically different states especially 
for small M. Both of these profiles are also clearly dis¬ 
tinguishable from that of an unscreened pair. The radius 
of a small screened pair with M = 2 is close to the ex¬ 
pected value 2^/2 + 7^ even for small system size N = 6 , 
which is significantly different from the radius R = 2 £ 
of the unscreened pair. We note that the smallest pair 
with M = 1 is not meaningful for the purposes of braid 
statistics, presumably due to the non-negligible overlap 
between the two quasiholes. Therefore, in our discussion 
below of the experimental feasibility in Sec. lYll we will 
only consider M > 2 . 
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V. MANY QUASIHOLES 


We have shown that the 2-quasihole wave functions 
can be produced in a system with two species of bosons 
interacting with a 3-body potential. Wave functions for 
states containing 2 n quasiholes, in principle, can also be 
produced as low/zero energy states of a system of N z- 
bosons and 2n auxiliary bosons under a 3-body inter¬ 
action Hamiltonian. Specifically, one can see from the 
bipartite form of the MR 2n-quasihole states (shown in 
Eq. (A 6 )) that they have zero energy for the interac¬ 
tion V = Yz 7 .z + Vzz-q- As we describe in the following, 
the zero-energy states of this Hamilonian have degenera¬ 
cies that exceed the 2"’“^-fold degeneracy of the FQHE 
quasihole states, due to the additional degree of freedom 
associated with the wave function of the auxiliary bosons. 
However when these degrees of freedom are frozen - for 
example by fixing the locations of the auxiliary bosons to 
produce localized quasiholes, the expected 2 "“^ degener¬ 
acy is obtained. 


The uninteresting degeneracies arising from bosonic 
edge excitations of the quantum Hall system can be re¬ 
moved by constraining the edge with a momentum cut 
off. We diagonalize the interaction V = -|- for a 

system of A = 6 z-bosons and four auxiliary bosons, with 
a cut-off on single-particle angular momentum equal to 
^max = (corresponding to the standard MR 4-quasihole 
wave function 'p 4 qh in a form of Eq. ([ 6 )). For simplicity 
we study only those sectors that have the same total an¬ 
gular momentum as the MR quasihole state. Numerically 
we find that the energy spectrum has three zero-energy 
states. 


The counting can be explained by noting that the zero- 
energy states should have the form 'I'(z, 77 ) x 3 ( 77 ) where 
4 '(z,r 7 ) is a 2n-quasihole state of the FQH system and 
is annihilated by Vzzz- Taking { 77 ^} as dynamic coordi¬ 
nates for the auxiliary bosons, '£'(^, 77 ) is also annihilated 
by Vzzrj because of the binding of the quasihole to the 
bosons. The function 3 ( 77 ) is any possible wave function 
of the 2n auxiliary bosons. A basis for space of such zero- 
energy states can be obtained as 77 ) x h^(r]) where 

ipF,x are the quasihole basis states introduced in Eq. ([^. 
hfj,{r]) are wave functions of 2n bosons which can be in¬ 
dexed by a partition (non-decreasing ordered sequence) 
/i of length 2n and has a form hf^{x) = S[xi^X 2 ^ ...]. 
The total angular momentum of this state is 


]\J-2 ^ ^ jp 

Lf,\f = (34) 


fc=i 


1=1 


Those states that have the same total angular momentum 
^ -I- N{n — 1) as the MR state of 2n quasiholes satisfy 

the constraint + Sf=i — !)• Correct 

degeneracy can be estimated by enumerating the states 
that satisfy this constraint. For the above mentioned case 


TABLE I: Number of zero-energy states for 2 — 6 quasiholes 


number of 
quasiholes 2n 

^max 

degeneracy for 
bosons 

degeneracy for 
localized quasiholes 

2 

A- 1 

1 

1 

4 

A 

3 

2 

6 

A + 1 

11 

4 


of 277 = 4 quasiholes, the states are 

(F = 0,A = [],77=[0,0,0,0]) 
(F = 2, A =[0,1], 77 =[0,0,0,2]) 
(F = 2, A =[0,1], 77 =[0,0,1,1]) 


In order to produce localized quasiholes, we fix the values 
of 77 i’s, and then the degeneracy of the zero-energy states 
reduces to the expected number 2"“^ = 2. It can be 
explicitly checked that the space spanned by the 2 "“^ as 
obtained states is identical to the space of the 2"“^ MR 
quasihole states with same quasihole positions. 

The same calculation for six quasiholes is performed 
with L = A^/2 + 2N and Imux = N + 1. There are 11 
zero-energy states for 6 77 i-bosons immersed in A = 6 z- 
bosons, and the degeneracy reduces to 4 when fixing the 
r]i positions. The degeneracy of 11 can again be obtained 
by enumerating {F, A, 77 ) that satisfy the mentioned con¬ 
straint. 

We summarize the degeneracy results in Table |T] and 
conclude that the MR 2 77 -quasiholes states can be pro¬ 
duced by turning on a 3-body interaction V = Vzzz + Vzzt] 
in a system of 2n auxiliary bosons plus A z-bosons and 
destroying the auxiliary bosons at some particular posi¬ 
tions, leaving 2n quasiholes at those locations. 

For larger values of Zmaxj the number of zero energy 
states is greater due to edge excitations. We have not 
studied the counting of such states. 


VI. FEASIBILITY IN ULTRA-COLD ATOM 
SYSTEMS 


A. Preparing the model wave functions 


FQHE states can in principle exist in rapidly rotat¬ 
ing gases of ultra-cold bosonic atoms. By an adiabatic 
scheme implemented by Gemelke et one can reach 
the ground state at a given total angular momentum L 
for A bosons. We have discussed in the Sec. Hill that the 
2-quasihole qubit states can be produced with a 3-body 
Hamiltonian. We now study how a Pfafhan wave func¬ 
tion of the form 'I'2qhTM(7?i) 772 ) can be prepared. We 
will only consider M > 2 cases as we discussed above 
that M = 1 is too small to show meaningful information. 
The standard Hamitonian in the LLL with a 3-body in- 
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teraction (in the rotating frame) is 

N+2 

H = Vo'^ 5{n-fj)S{n-fk) + {uJz-^)L;^ + (u;r,-^)Lrj, 

i,j<.k 

( 35 ) 

where is the rotation frequency, Wz and ujri are the 
harmonic confinement frequencies, and and are 
total angular momentum operators for the z and r] parti¬ 
cles respectively. We assume the interaction is indepen¬ 
dent of the species which produces the 4'(i)(2) state as 
its unique zero-energy state at Lq = iV^/ 2 . As shown 
above, the other 2-quasihole state ^'(i)(2) — 2d>(i2)() can 
be produced as the first excited state for an appropriately 
chosen interaction. For L = Lq + M, 'I'(i)(2)-Pm(??i,?? 2) 
has zero interaction energy, but it is in general not the 
only zero-energy state; other zero-energy states can be 
constructed wherein the additional angular momentum 
M is absorbed by z particles. We ask if the edge ex¬ 
citations of the z particles can be suppressed by tak¬ 
ing ujz > ojjj, following the strategy used in Ref. UHl 
The exact ground state of this Hamiltonian at M = 2 
and M = 3 is well approximated by ^'(i)(2)^’m(?7i, ^72), 
while for M > 4 the second or third excited state has 
a pair correlation function matching well with that of 
'I'(i)(2)^’m(??i, »72)- At minimum, the state of M = 2 and 
M = 3 can be adiabatically prepared in principle. Then 
the measurement for pair correlation function can be im¬ 
plemented through single-atom detection combined with 
a short time-of-flight expansion. More detailed discussion 
can be found in Ref. [901 


B. 2-body contact interaction 

We have shown above that many relevant topological 
states with two or more quasiholes can be engineered with 
3 -body interaction. However, the 3 -body interaction is 
not easy to implement. A more realistic interaction in 
cold atom systems is the 2-body contact interaction: 

H = Vo^ 5 {ri-rj). ( 36 ) 

i<j 

The MR Pfafhan state is not the exact ground state of 
this Ha milton ian. Encouragingly, however, it has been 
showiP^tSSEIl that this interaction produces an incom¬ 
pressible state at = 1, which is accurately, though 
not exactly, described by the Pfafhan wave function. We 
therefore ask to what extent the above results carry over 
to the 2 -body contact interaction. We focus on the fol¬ 
lowing two aspects: (i) Can a 2 -body interaction pro¬ 
duce the low-energy Pfafhan quasihole states with the 
expected degeneracy? (ii) In order to detect statis¬ 
tics through pair correlation function measurement, we 
ask whether a 2-body interaction can produce a ground 
state of which the pair correlation function for two aux¬ 
iliary bosons is similar to that of the Pfafhan wave func¬ 
tion 4>2qhC’M(?7i,f?2)- We will only discuss the 'I'(i)(2) 


state for simplicity because similar conclusions apply to 
'f'(i)(2) — 24 '(i 2 )(). We will also focus on M > 2 cases as 
M — 1 has been shown ineffective in rehecting statistics 
information correctly. 

To answer the hrst question, we perform exact diag- 
onalization of 6-10 bosons with an additional 2 and 4 
auxiliary bosons within a Hilbert space specihed by a to¬ 
tal angular momentum L and a cut-off on single-particle 
angular momentum Zmax- The value of L is hxed accord¬ 
ing to the Pfafhan quasihole wave functions, for example, 
L = N ^/2 for 2 -quasihole states 'I'(i)(2) and 'I'(i2)(). The 
value of Zniax is taken to be equal or greater than the 
largest single-particle angular momentum in wave func¬ 
tions. In order to tune the 2 -body interaction, we write 
the interaction potential as 


V — Vzz + Ait 4 „ -I- A2I4, 


( 37 ) 


where Ai and A2 control the relative strength of z-rj 
and 77-77 interactions. Because the 3 -body interaction 
V = Vzzz + VzzT] excluding 77-77 interaction term {Vztjti) 
produces degenerate Pfafhan quasihole states, we keep 
A2 = 0.0 and only tune the value of Ai for the 2-body in¬ 
teraction. We have not been able to hnd any values of Ai 
and Imax that produce a degenerate or quasi-degenerate 
set of low-energy states with signihcant overlap with the 
Pfafhan quasihole wave functions, for either two or four 
quasiholes. 

Now we study the second question for M = 2. We 
recall that the Pfafhan wave function {i){2)PMi'nijV2) 
can be approximately produced as the ground state of the 
Hamiltonian in Eq. ( 351 with > w,, at L = /2 + M. 

In fact, if we take /max to be the minimal value /max = 
N — 1 , the zero-interaction-energy states have the ex¬ 
act unique form 'i 'without requiring the 
conhnement term in the Hamiltonian. We now diagonal¬ 
ize an isotropic 2 -body interaction (Ai = A2 = 1 . 0 ) with 
the same /max = — 1 and L = A ^^/2 -|- M for M = 2 . 

The pair correlation functions for the 77-bosons of the 
ground state are clearly different from that of the wave 
function '^{i){2)PM{vii'n2), as shown in Figure|^(a). We 
next tune the parameters Ai_2 and hnd that A2 = 10 and 
Ai = 1.6 give good agreement for system size of A^ = 8, 
which also apply to A^ = 6 and A^ = 10 as shown in 
Figure]^ (b). 

For M > 2 , a conhnement potential that is slightly 
stronger for z-bosons than 77-bosons needs to be included 
in the Hamiltonian in order to help the rji ’s absorb all of 
the additional angular momentum to produce the multi¬ 
plicative factor Pm ivij''12)■ We denote the difference in 
the harmonic conhnement frequencies of the two species 
of bosons as Sui = Wz — By diagonalizing the new 
Hamiltonian with different at L = / 2 -\- M and 

/max = A^ — 1, we hnd that the pair correlation func¬ 
tion of the ground state depends sensitively on the value 
of 6 u}. Figure (c) shows the pair correlation function 
of the ground state obtained with parameters Ai = 1.6, 
A2 = 10 and 5 uj = 0.15 at M = 3 for different N. It is 
not possible to hnd a single set of parameters that gives 
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FIG. 6 : (Color online). The pair correlation functions g{r,0) 
of the 77 -bosons of the ground state obtained from diagonaliz¬ 
ing the 2-body Hamiltonian with parameters shown in the fig- 
nre are plotted for different system sizes N = 6 (blue square), 
N = 8 (green circle), = 10 (red diamond). As comparison, 
the pair correlation functions calculated from the model wave 
functions ^'{i){ 2 )F’m( 7 ?i, 772 ) are also shown with solid lines. 
The plots for N = 8 and 10 have been shifted up by 0.05 and 
0.1 units respectively for ease of depiction. (a)(b) For M = 2, 
an isotropic contact interaction fails to produce a similar pair 
correlation function as that of the model wave fnnction, but 
it is possible to find a universal set of interaction parameters 
Ai ,2 for different system sizes that gives agreement, (c) For 
M = 3, a slightly different confinement potential for z- and 
77 -bosons is necessary to produce a desired ground state. How¬ 
ever, the proper value of parameter Slj sensitively depends on 
the system size. 


the desired pair correlation function for all N. 

Up to now we have studied the state; the same 

arguments apply to the other state ~ 2 'I'(i 2 )() with 


optimal parameters Ai = 0.4, A 2 = 10 and 6u! = 0.023 
for N = 8. In the study of M = 2,3 cases, all exact 
diagonalizations are performed with the smallest value 
of ^max as specified by the wave function. A larger /max 
generally requires a different set of parameters. In con¬ 
clusion, 2-body interaction fails to produce a realization 
of the Pfafhan quasihole wave functions with the desired 
degeneracy, at least for all the systems we have consid¬ 
ered. 

VII. CONCLUSIONS 


We have pursued the idea that the non-Abelian quasi¬ 
holes in the Pfafhan fractional quantum Hall states can 
be captured and manipulated by introducing auxiliary 
bosons of a different species. We have shown how, as¬ 
suming a precise control of the Hamiltonian, this can 
be used to produce qubit states with desired topological 
properties. Furthermore, a measurement of the pair dis¬ 
tribution function of the auxiliary bosons can provide a 
direct conhrmation of the Abelian and non-Abelian braid 
statistics of these quasiholes. 
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Appendix A: Equivalence of Bipartite and Pfaffian 
Representations of 2n-quasihole states 

We demonstrate that the Pfaffian wave function of 2n- 
quasihole states can be represented as an equivalent bi¬ 
partite form. We hrst consider the 2-quasihole states for 
simplicity. There are two 2-quasihole states in the bi¬ 
partite representation: one in which the two quasiholes 
are constructed in the same partition and another one 
in which the two quasiholes are constructed in different 
partitions. They have the form: 


^(12)()(^,ffi,??2) = 5 < W{z2i - ■ni){z2i - m) n - Z2:j-lf > , 


^<j=l 


(Al) 


(l)( 2 ){z,r]i,r] 2 ) = 5 <1 W{Z 2 ^ - ?7l)(Z27-l - V 2 ) - Z 2 Jfiz 2 ^-l - Z 2 j-lf 


2=1 




(A2) 
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Following, we demonstrate that these states can be repre¬ 
sented in terms of the Pfaffian of a matrix. We will take 
'k(i 2 )() as an example and the proof for can be 

completed by simply replacing Z 2 i in the factor {z 2 i — r] 2 ) 


with Z 2 i-i throughout the derivation. 


Factorizing out a Jastrow factor J(z) = Y\^^j{z^ — Zj) 
from the wave function gives 


^ii2)oiz,VuV2) = J{z) X W(.Z2i - rii){z2i - V2) 


= (-I)"! 2 1) J(z) X ^ W{z2i - r]i){z2i - ?72)det 


W<j=l{z2i - Z2j){z2i-1 - Z2j-l) 

N 

Wj = liz2i - Z2j-l) 

1 


2=1 


^2i — ^2j-l 


(A3) 


= (- 1 ) 4(2 i)j(2)x sgTL{a)Al\\{z2^-m){z2^-r]2)X\—^ 




N/2 


Z2t - Z2cri-1 


The symbol A denotes antisymmetrizatio n ov er all coor- has been used to obtain the second line of Eq. Ad The 
dinates Zi, and an identity due to Cauchj®^, third line comes from the expansion of the determinant 

where the sum is computed over all permutations a of 
1 Z 2 j){z 2 i-i -Z 2 j-i)t^e set {l,2,...,iV/2}. Note that all the permutations 

det--- = (—1) 4 '■ 2 •'- - -yield the same contribution, thus the wave function can 


- Z2J-1 


niy(^ 2 i — Z 2 j-i) be simplified as 

(A4) 


'^ii 2 )o{z,VuV 2 ) = (- 1 )”^" ^\n/ 2 )\J{z) X 


i=l 


= (-l)f (f-i)(7V/2)!J(z) X Pf 


{z2i - m)iz2i - m) 

Z2i — Z2i-1 

{zi - r]i)izi - m) + (i ^ j) 


Zi - Zi 


(A5) 


which is identical to the Pfaffian representation in Eq. 
(10 1 up to an overall normalization factor. The defini¬ 


tion of Pfaffian has been used in the last step, and the 
Gaussian factors are suppressed for simplicity. 

The above derivation can be easily generalized to 

K 

2 n-quasihole state by replacing the factor Y\i^i{z 2 i — 


■ni){z2i-i - m) in Eq. (A2) with niLi[(^ 2 i - Va){z 2 z - 


r]p)...{z 2 i-i - r]p){z 2 i-i - Vcr)---]- Thus a bipartite rep¬ 
resentation for the Pfaffian wave function in Eq. (|^ is 
constructed: 




(a/3...)(p<T...) = 5 < W{Z2^ - ria){z2r - r||i)...{z2^-l - r]p){z2i-l - ?7^)... (Z2i “ ^2j)^(z24-l “ Z2j-lf 


(A 6 ) 


i=l 


i<j=l 


Appendix B: 2-quasihole representation of the 
4-quasihole qubit states in certain limit 


states in Eq. in the limit that two of the quasiholes 
are sent to infinity. We take = r, 774 = re*® and r —>■ 


We now show that the 2-quasihole states 'I'(i)( 2 ) and 
'f'(i)( 2 ) ~ 24'(i 2)() can represent the two 4-quasihole qubit 
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oo for Eq. ((^. Now 'I'(i 3 )( 24 ) and d>(i 4 )( 23 ) reduce to 
the standard MR 2-quasihole wave function 4 '(i)( 2 ) with 
quasiholes located at rji and r] 2 , while 4 '(i 2 )( 34 ) reduce to 
'£'( 12 ) 0 , with a same normalization factor for all of them 
as shown below: 

^(13)(24),^'(14)(23) 

T N 

'f'(12)(34) r 6 2 W(i 2 )()- 

In this limit, the factor a; —>■ 1 -|- (1 — ) = 

1 + e(ryi — 772 ), where e = ^— is a very small num¬ 
ber. Before we apply any approximation, we rewrite 


^(13)(24) ± V^^'(14)(23) in Eq. § as 

{x±^J 2;)4'(14)(23) + (1 — a^)'I'(12)(34) (B2) 

with the help of the identity 

'i'(12)(34) ~ ^(14)(23) = “ ('i'(12)(34) ~ 'i'(13)(24)) (B3) 

from Ref. |25l By applying ^/x —>■ l-b |e(?7i — 772 ), the two 
topologically different 4-quasiholes wave functions v[/(i.V’) 
in Eq. ([^ are written as (for ttt = 1): 


(1 ± (1-b |e( 77 i - 772 ))”^^'^ v^)^(i 4)(23) + (1 - a^)^(i2)(34))- (B4) 


Thus 


^ {^rji - 772 )^/®4'(14)(23) (m - 

and 

(•_ ~ ^2)(4'(i 4)(23) “ 24'(i2)(34) )] (? 7 l “ ^2)®'^® (4'(i) (2) - 24'(42)()) 


(B5) 


(B6) 


Therefore reduces to 4 '(i)( 2 ), and reduces to a 
linear combination of the two kinds of 2 -quasihole wave 
functions ( 4 '(i)( 2 ) — 24'(i2)()), with prefactors (771 — 772 )^/® 
and (771 — 772 )®^® respectively. 

In the above derivation we have neglected the factor 


_ y->4 |t7q I 

e jT;q_ ( 5 )^ This factor simply contributes 

to an overall normalization constant of the wave func¬ 
tion and is only considered when the dependence of the 
normalization constant on the locations of quasiholes is 


discussed, as in the derivation of Eq. (30) in Sec. IV B 


In that case, a more precise version of Eq. (B51 and (B 6 ) 
is needed as 


4/(5 = (rh 


??2)^^®4'(i)(2)e 


1^1 \‘^ + \V2 P 


772)®/®(4'(i)( 2) - 24'(i2)())e 


hi l^ + h 21 ^ 

sP 

(B7) 
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